Universality versus material dependence of fluctuation forces between metallic wires 
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We calculate the Casimir interaction between two parallel wires and between a wire and a metall 
plate. The dielectric properties of the objects are described by the plasma, Drude and perfect 
metal models. We find that at asymptotically large separation interactions involving plasma wires 
and/or plates are independent of the material properties, but depend on the dc conductivity a for 
Drude wires. Counterintuitively, at intermediate separations the interaction involving Drude wires 
can become independent of a. At smaller separations, we compute the interaction numerically and 
observe an approach to the proximity approximation. 



Effective interactions between nanowires and nan- 
otubes have attracted lots of attention due to their grow- 
ing applications in micro- and nanomechanical systems 
[IHl]. The knowledge of the interactions between single 
walled carbon nanotubes (SWCNT) with different chi- 
rality and hence electromagnetic response is important 
to separate a polydisperse solution of SWCNT in frac- 
tions of equal chirality [1]. Under many circumstances, 
van der Waals or Casimir forces are the dominant in- 
teraction and hence a precise understanding of them is 
needed. Furthermore, cylindrical shapes are important 
for precision Casimir force measurements, in comparison 
to spheres, because of the larger effective area of interac- 
tion [5J |B] . Approximations of the Casimir force between 
cylinders and plates [5] have shown that the tempera- 
ture dependence varies based on the description of the 
material properties. Thus there is a need for exact calcu- 
lations of the Casimir force for cylindrical shapes taking 
into account realistic material response. 

It has been demonstrated that Casimir interactions 
strongly depend on the combined effects of shape and 
material properties, see, e.g., [3 E]. The interplay is 
particularly strong for quasi one-dimensional conduct- 
ing materials due to strongly anisotropic collective charge 
fluctuations. Indeed, for two parallel perfectly conduct- 
ing wires of distance d the retarded interaction energy per 
length is £/L ^ hc/cP, apart from a logarithmic factor 
[9] . It decays only slowly compared to the retarded inter- 
action £ / L ^ hcR^/d^ between insulating cylinders that 
do not support collective fluctuations. Most studies of in- 
teractions between one-dimensional systems over a wide 
range of separations concentrate on these two situations. 
However, low dimensionality in combination with flnite 
conductivity and plasmon excitations should give rise to 
interesting new effects that might be probed experimen- 
tally using, e.g., the coupling to mechanical oscillation 
modes. The often employed technique for these effects, 
the proximity force approximation (PFA) cannot capture 
the correlations of shape and material response since it is 
based on the interaction between planar surfaces. There 
have been attempts to compute the van der Waals in- 



teraction between cylinders (and plates) for particular 
frequency dependent permittivities [5l415j. However, the 
interplay between shape and material effects is not trans- 
parent in these works as they are limited either to perfect 
metals or to asymptotic limits. 

In this Letter, we employ the scattering approach to 
study the interaction between two infinitely long, parallel 
metallic wires and a wire and a plate that are described 
either by the plasma or the Drude dielectric function. 
We model the wires as circular cylinders of radius R and 
obtain analytical results for the interactions at distances 
much larger than R. We find regimes of different am- 
plitude and power-law, depending on the relation of the 
radius to the length Ao- = 27rc/(T with conductivity a and 
the plasma wave length Ap. 

Most interestingly, we find that the interaction involv- 
ing Drude wires approaches the universal interaction be- 
tween perfect metal wires at intermediate separations 
and becomes non-universal (material dependent) at dis- 
tances d > R^/Xcr- This behavior is explained in terms 
of the size of collective charge fluctuations in a Drude 
metal. For wires that support plasma oscillations, the in- 
teraction does show universality at asymptotically large 
distances. In order to confirm the validity of our large dis- 
tance expansions and to compare with the PFA at short 
distances we have performed numerical computations of 
the interaction. 

To investigate the Casimir interaction between a cylin- 
der parallel to a plate or another cylinder, we employ 
the scattering approach, which allows us to calculate the 
energy over a large range of separations between the ob- 
jects. Because of the translational symmetry along the 
cylinder axis, the Casimir interaction can be written as 
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with L the length of the cylinder, n the Wick rotated 
frequency, the wave number along the cylinder axis, 
and 1 the identity matrix. The matrix N factorizes into 
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the object's scattering amplitudes (T-matrices) that en- 
code the material dependence and distance dependent 
translation matrices that described the coupling between 
multipoles on distinct objects. 

Parallel wires - For two infinitely long parallel cylin- 
ders with radii R, aligned along the z-axis, the elements 
of the matrix N for polarizations a, /3 — E, M and partial 



waves m, m are 
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with T the T-matrix of the cylinder (see below). The 
translation matrix, U^'^ , relates the regulare vector waves 
M^;/„ = {n^ + kl)-^l'-V X [/„(pp)e*("«+'=^-)] and N^,f„ = 
K^^V X M^''^ with p = + k1 in the cylindrical co- 
ordinates (p, 9, z) of one cylinder to the outgoing vector 
waves of the other cylinder that are given through re- 
placing the Bessel functions /„ by the Bessel functions of 
second kind, if„. The elements are given by 



(3) 



where d is the distance between the centers of cylinders. 
Note that the translation matrix conserves polarization 
and thus is independent of the polarization index. The T- 
matrix for a cylinder with dielectric function e{icK) and 
magnetic permeability ii{icK) is diagonal in the partial 
wave number n but couples polarizations. Its elements 
are 
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wi th P„ = {nKJ{^R^K)){l/p'^ -l/p"), p' = 
^ e^K? + k1 and 
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A2,n is obtained from Eq. ([6|) by interchanging e with /i, 
and Aa^n follows from Eq. TSl) by replacing K'^ with 
and with J„. The elements Tf^^ are given by T^^ 
after interchanging e with fj,. Finally, antisymmetry in 
polarization yields T,f ^ = -T^*^. 

In order to investigate the impact of the material prop- 
erties on the Casimir interaction, we consider plasma, 
Drude and perfect metal cylinders with the magnetic per- 
meability fi = 1. The Drude dielectric function is 



e(ic/t) = 1 + 



(ApK)2 + A,k/2 
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and reproduces the plasma model for — 0. Using 
Eq. ([7| , the asymptotic behavior of the T-matrix element 



of Eq. Q for n = at small frequencies (k <^ 1, kz/K 
fixed) is given by 
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with C{k) « X/k^ /(27r2) if i? < Xp, C{k) = X„k/{4tt^) 
if K <C Xcr/Xp, l/Xa, and C{k) = for plasma, Drude 
and perfect metal cylinders, respectively. At small fre- 
quencies K but fixed kz/K, one has for Drude cylinders 
~ K, while for plasma and perfect metal cylinders 
one has T^^ - 1. Since T^^^' - k^, T^^'^ = T^'^ = 
and higher order elements associated with n ^ scale as 
K^l"!, it is sufficient for the large distance interaction to 
consider only the element T^^ ■ 

The Casimir interaction between metallic wires is in 
general complicated and no simple analytical expres- 
sion that applies to all length scales can be obtained. 
However, using Eqs. ^ and ^ along with T^^ given 
in Eq. ([s]), the asymptotic interaction at large sepa- 
rations, d ^ R, can be calculated in various limiting 
cases. For metals with diverging response at zero fre- 
quency, one usually observes universal behavior for the 
interaction at large separations. Indeed, that is what 
we obtain for the interaction between plasma wires (or 
plasma and perfect metal wires) which then becomes 
£ = -~hcL/{8TTd'^\n^id/R)) [HIIS]- This universal form is 
only applicable beyond an exponentially large crossover 
length d ~ Rexp{Xp/ R"^). Below this scale the inter- 
action becomes material dependent. For two plasma 
wires or a plasma wire and a perfect metal wire, with 
Xp/R > 1, the energy scales as -R/ {Xpd^ Iv?^"^ {d/ R)) . 
For numerical coefficients, see Fig. [l|a). 

For configurations involving at least one Drude cylin- 
der we find a rather distinct behavior that deviates from 
naive expectations for universality. For large distances 
d ^ R, Xcr we obtain two different scaling regimes that 
are separated, up to logarithmic corrections, by the curve 
d/R ^ ^/d/Acr, see Fig. [l|b). The unexpected feature 
is that the interaction is universal in the regime where 
d <^ R^/Xa- If the distance is increased beyond this 
crossover scale (with all other length scales kept fixed, 
see arrow (1) in Fig. [T]Jb)), the interaction becomes ma- 
terial dependent and up to logarithmic corrections, scales 
as —R?/{Xryd^) for a Drude wire interacting with another 
Drude wire or a plasma or perfect metal wire. For de- 
tailed forms of the interactions in this limit, see Fig.[l|^b). 
If, however, the radii of the wires are increased in the 
same way as their distance {d/R fixed, see arrow (2) 
in Fig. [ij^b)), finite conductivity becomes unimportant 
at large distances and the interaction assumes the per- 
fect metal form. An intuitive explanation of this non- 
universal large distance behavior is given below. Note 
that the decay of interactions between insulating wires 
for c? ^ i? scales as hcLR'^/d^ with a material dependent 
coefficient. 
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FIG. 1. Summary of the different forms of interaction be- 
tween two wires and a wire and a plate. Shown are the 
rescaled interaction energies per cylinder length, £/{hcL). (a) 
Interaction involving a plasma (P) wire with another plasma 
wire, a perfect metal (PM) wire or a plate. The asymptotic 
results apply sufficiently far away from the separating curve 
ln(d/i?) ~ {Xp/Rf and for d/R, Xp/R > 1. (b) Interaction 
involving a Drude (D) wire with another Drude wire, a plasma 
wire, a perfect metal wire or a plate. The separating curve is 
given, up to logarithmic corrections, by d/R ~ ^J~df\,. The 
shown expressions hold for d/_R, d/Ao- 2> 1 and d ^ Ap/Ao-. 
Depending on the relative size of length scales, different 
regimes can be reached: Dashed arrow (1) corresponds to 
an increasing distance d which ultimately leads to a non- 
universal interaction. Dashed arrow (2) indicates an overall 
increase of the geometry (i.e., djR fixed) with constant con- 
ductivity leading to a universal interaction. 



Wire parallel to a plane - Now we consider a wire that 
is parallel to a plate. We assume that the plate is in the 
y — z plane and its distance from the cylinder (2) axis is 
d. Then the matrix N in Ec^. ([T]) can be written as 



n: 



. „ .. . -■. 



in'=E,M 



2 + k2 



(9) 



where ky is j/-component of the wave vector, kj^ 



{ky^kz), and the matrix Dmk^'y,\!.±Y converts between 
plane and cylindrical waves |16j . Further, T^^ are the 
T-matrix elements for the plate which are given by the 
usual Fresnel coefficients [TB]. For perfect metal plates 
Tip, — T^^ — 1 and for small k at fixed k±/K one has 



for the plasma model Tjf^ 
the Drude model 



k_L ~ 

^ -^M=l + 0{XpK) and for 
l + 0{XaK). Due to this 
identical behavior of the plate's T-matrix at small k, the 
interaction at asymptotically large distances is indepen- 
dent of the model that describes the metal plate. 

We calculate the Casimir energy between a wire and 
a plate at large separations (d/R ^ 1), using Eqs. ([T]), 
^ and ([9|. We again find two different scaling regimes 
that are separated by curves that are given by the same 
expressions that we found for two wires, see Fig. [l] In 
the universal regime (perfect metal), the interaction is 
£ = ~hcL/il6Trd'^\ii'^{d/R)) US] for a plasma wire 
at asymptotically large d ^ R and for a Drude wire at 
intermediate distances with A^, Ap/Ac <C d ^ R'^/Xa, 
see Fig. [IJb). In the other regime the interaction is non- 
universal with the energies up to logarithmic accuracy 
scaling as f ^ hcLR/ (Xpd^) for a plasma wire and £ ~ 
hcLR? / [Xad^) for a Drude wire. For the precise form see 
Fig. [l|b). This can be compared to the faster decay that 
is observed for an insulating wire interacting with a plane 
which has £ - hcLR^/d^ [l6]. 
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FIG. 2. Ratio of the numerically computed force and the 
universal (squares) and non-universal (circles) asymptotic re- 
sults of Fig. [l] for two Drude wires with Xp/R = 0.5 and 
Act = Ap/27.4. The inset shows the ratio of the numerical 
force to the proximity force approximation (PFA) for Drude 
cylinders with Xp/R — 0.5 (diamonds) and Xp/R — 0.05 (tri- 
angles). 

Numerical results - It is interesting to compare the 
asymptotic results with exact numerical calculations to 
identify the regions in which the asymptotics are correct. 
We calculate numerically the integrals and determinant 
in Eq. ([T]). The infinite matrix N given by Eq. ^ is 
truncated at a finite number n = rXmax of partial waves. 
JT-max depends on the separation between the objects and 
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is chosen such that the energy varies by less than 0.01% 
upon increasing Umax by 10. As the separation becomes 
shorter, rimax increases. While for d/R = 10, rimax = 9 
is sufficient, we need to use nmax = 91forc?/i? = 2.1. 

Figure [2] shows the ratio of the numerically computed 
force between two Drude wires and the corresponding 
asymptotic results (universal and non-universal regimes, 
see Fig. [ijb)) versus ln(2c?/_R). The material parame- 
ters are chosen as Ap/i? — 0.5 and Ao- — Ap/27.4 which 
correspond to gold with Ap = 137 nm and A^- ~ 5 nm. 
At intermediate separations, the force normalized to the 
universal result approaches unity whereas at asymptoti- 
cally large separations the force normalized to the non- 
universal result tends to unity. This confirms the validity 
of the crossover shown in Fig.[ljb). We also compare our 
numerical results for the force at short separations with 
those obtained from the proximity force approximation 
(PFA) which is based on the Lifshitz formula for parallel 
plates made of the same material as the wires. The in- 
set in Fig. ([2]) shows the ratio of our numerical result for 
the force and the PFA result versus d/R for Drude wires 
with Xp/R — 0.5 and Ap/i? = 0.05. Our data support the 
validity of the PFA in the limit of vanishing separations. 

In summary, we have calculated the Casimir force be- 
tween two metallic cylinders and a metallic cylinder par- 
allel to a plate. We find that a significant feature of 
the interaction between a Drude Cylinder with another 
Drude cylinder or a plate is that upon increasing the sep- 
aration, the interaction can move from a universal regime 
to a non-universal one. This behavior can be under- 
stood from the wave equation for the electric field inside a 
Drude cylinder. For imaginary frequencies u = icn, the 
Helmholtz operator + e(w)(a;/c)^ for a good Drude 
conductor becomes — Stt^k/Ao-. We are interested in 
the maximal scale of the field and hence charge fluctu- 
ations for a given k. With the smallest transverse wave 
vector kx, ky ^ 2tt / R we find the dispersion relation 



\k,\^ R~^^k/k,~1, k.^XJR^ 
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Hence, collective charge fiuctuations on arbitrarily large 
scales exist only for k > Kc which is a consequence 
of dimensionality that does not appear in the absence 
of transverse constraints {R — )■ 00). For k < Kc 
charge fluctuations break up into clusters of typical size 
~ R/^/l — k/kc due to flnite conductivity. The spec- 
tral contribution to the interaction between cylinders at 
distance d is peaked around k ^ 1/d. If d < 1/^^c 
{d/R < ^d/\fj, see Fig. [l|b)), collective charge fiuc- 
tuations contribute strongly to the interaction and ren- 
der it universal as for perfect metal cylinders for which 
Kc ^ 1/cr — > 0. In the asymptotic regime with d > 1/kc 
{d/R > i/d/ACT, see Fig. [l|b)), finite conductivity pre- 
vents fluctuations on arbitrarily large scales and hence 
the interaction is proportional to cr, i.e., non-universal. 
It is important to note that as R goes to zero, Kc be- 



comes larger, and in consequence the flnite conductivity 
of cylinder becomes more important. 

Finally, an estimate of the interaction between two 
gold wires based on the asymptotic expressions in Fig. [l] 
with R = lOnm, length L — 100/im, Ap = 137nm and 
Act — 5nm at a distance d = 200nm yields a force of 
ss IpN within the plasma description and « 27pN within 
the Drude model. These forces are experimentally de- 
tectable and allow a much clearer distinction between 
plasma and Drude model predictions as compared to two 
plates or a plate and sphere [S]. This is of particular 
importance in view of recent experimental flndings that 
interactions between metals might not be consistent with 
the Drude model [T7] . 
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